On resumming periodic orbits in the spectra of integrable systems 
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Spectral determinants have proven to be valuable tools for resumming the periodic orbits in the 
Gutzwiller trace formula of chaotic systems. We investigate these tools in the context of integrable 
systems to which these techniques have not been previously applied. Our specific model is a strobo- 
scopic map of an integrable Hamiltonian system with quadratic action dependence, for which each 
stage of the semiclassical approximation can be controlled. It is found that large errors occur in the 
semiclassical traces due to edge corrections which may be neglected if the eigenvalues are obtained 
by Fourier transformation over the long time dynamics. However, these errors cause serious harm 
to the spectral approximations of an integrable system obtained via the spectral determinants. The 
symmetry property of the spectral determinant does not generally alleviate the error, since it some- 
times sheds a pair of eigenvalues from the unit circle. By taking into account the leading order 
asymptotics of the edge corrections, the spectral determinant method makes a significant recovery. 

PACS numbers: 03.65.Sq, 02.30.Ik, 02.30.Lt 



I. INTRODUCTION 



It has been generally accepted that spectral determinants, or zeta- functions, provide optimal semiclassical estimates 
of the individual energy levels of classically chaotic systems [0, ^, || . These resummations of the periodic orbits in 
the Gutzwiller trace formula |^, originally developed for time-independent systems, may be obtained in a variety 
of ways. The formulation of Berry and Keating |^ relies explicitly on the instability of the periodic orbits in such a 
way as to give an expression that does not depend on a sharp cut-off of the orbit periods. In contrast, Bogomolny's 
approach |^ ||] reduces the problem to the quantization of a map over a Poincaree surface of section, without making 
any assumption about the nature of the classical motion. However, two sharp boundaries are introduced. On the 
one hand, the Poincaree section itself is bounded if the constant energy surface is compact. The limited area of the 
section, which corresponds to a finite dimension of the Hilbert space, then leads to a sharp cut-off in the period of 
the orbits. The recent paper of Eckhart and Smilansky |^ also works with a quantum map in a bounded region, but 
this is obtained stroboscopically instead of with a Poincare section. 

The classical motion of integrable systems is restricted to invariant tori, which determine closed invariant curves 
(also tori) of the Poincare mapping. For integrable systems. Berry and Tabor jioj established the equivalence of the 
Gutzwiller trace formula to the general forms of Bohr-Sommerfeld quantization. The latter method is evidently the 
most efficient for the calculation of individual levels. Nevertheless the exercise of showing the equivalence of both 
methods had the merit of clarifying the role of periodic orbits (forming continuous tori for higher dimensional systems) 
in the density of states. Indeed it raises an important question: The Berry- Tabor equivalence involves the complete 
set of periodic orbits, so how might one obtain correct energy levels from the resummation of a finite selection of short 
orbits? 

Our interest in this paper is resummation approximations for the spectra of simple integrable systems. Particularly, 
we will focus on the effect of introducing boundaries that do not interfere with the local classical tori. For simplicity, 
consider such a system with a single degree of freedom. For the classical Hamiltonian H{I), where (1,9) are the 
action-angle variables ||ll|, the Bohr-Sommerfeld levels are 

En = H(^h(^n+^^^ . (1) 

Actually, it suffices to consider a case for which the action variables are those of the harmonic oscillator 

I^lip'+q'), (2) 

rendering the semiclassical approximation of (|^) exact. It is important to point out that large-n states are well 
approximated by (|l|) even for a general nonquadratic dependence of the Hamiltonian on the phase space variables. 
The main question is how well a spectral determinant method converges to the eigenvalues E„ in this case. 
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A conceptually clean approach to answering this question follows in the spirit of Creagh [|2|, who applied an 
extension of resummation methods |p^ , ^ to the perturbed cat maps. He found worsening errors in the semiclassical 
trace formula as the perturbation brought the map out of the hyperbolic regime, but was unable to continue the 
perturbation all the way to an integrable regime. Here, since / is a constant of the classical evolution, our consideration 
can be limited to a ring I_ < I < X+ , thus obtaining a finite phase space and hence a finite Hilbert space of dimension 
N = {(X-f— where { • } denotes the integer part. A stroboscopic map of period r can be defined and quantized 
in such a way that the dynamics is determined by the evolution operator 

un= ^''pi^J^"'") ■ 

n=N- ^ ^ 

The quantity hN^ (^A^+) is the smallest (largest) quantized action greater (smaller) than I_ (X+) and \n) are the 
eigenstates. The definition of the spectral determinant as det (1 — zUn) gives an A^*'* order polynomial in the variable 
z and the roots of the characteristic equation 

PAr(z) = det (1 - zUn) = (4) 

are just 

Zn = exp (^^EnT^ . (5) 

Thus, knowledge of the phases 4>n of the zeros of the spectral determinant, or zeta-function, determines the eigenvalues 
of the stroboscopic map within the ring. These coincide with a subset of the eigenvalues of the original system 
possessing an infinite Hilbert space. 

The eigenenergies of the original system can be inferred from 

= -(0n + 2fc7r) . (6) 

r 

The determination is unique if, for the energy interval of the ring (-E-, i?+) — (i?(X_), E(T+)), the constraint 

{E+-E^)^<27T (7) 

is respected. This scheme for locating eigenvalues of the original Hamiltonian resembles that of but our objective 
of comparing semiclassical methods can be achieved also in the context of the continuous dynamical system confined 
to the ring, or even the discrete map defined by (|^). 

In the semiclassical limit, Tr(J7^) is expressed as a sum over periodic orbits of period I for the classical stroboscopic 
map generated by the Hamiltonian H(I) in the time r within the ring. These orbits, with action such that 

dH(I) 27rm 

define continuous curves, unless 7 = 0. To simplify the treatment, we consider X_ > 0, which excludes the isolated 
periodic orbit at the origin. 

Note that the periodic orbits of the stroboscopic map comprise only a small subset of those of the original system 
with continuous time. Indeed, all the orbits of the latter are periodic whereas only those full orbits with periods that 
are rationally related to the stroboscopic time, r, are periodic in the map. The map resembles a Poincare map of an 
integrable system with two degrees of freedom in the sense that all orbits lie on invariant curves, but the curves made 
up of periodic orbits form a dense set of zero measure. 

The semiclassical form for the spectral determinant is obtained from the expression 



det (1 — zU at) = exp 



1=1 



(9) 



Even though the series only converges for |z| < 1, it is possible to follow [jisj in noting that the Taylor expansion of 
the exponential in Eq. (|p) can be identified with the finite expansion of the spectral determinant (^): 

PNiz) = 1 + CIZ + C2Z'^...+CNZ'^ , (10) 
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where the coefficients are given by the following recurrence relation 

i^c,_,Tr(C/^). (11) 



Ck 



k 



In this way the periodic orbits up to period I ~ N determine, in principle, all the eigenvalues of Un- 

A further reduction of the period of the orbits used in the semiclassical spectral determinant results from the 
symmetry relation for the coefficients [Q 

Ck = dct{UN)cN^k . (12) 

The coefficients containing traces for long iteration times, I > N/2, can be obtained from the coefficients with I < N/2. 
Hence, only the orbits of period / < N/2 are needed. This is a particularly simple example of "bootstrapping" |l^ 
resulting from the finiteness of the Hilbert space. For chaotic systems, the symmetry ( p^ is a fundamental tool 
because the number of periodic orbits increases exponentially with I. This is not so for integrable systems. 

In Section ||, we derive the semiclassical formula for Tr(J7') of the integrable map defined in (||). It will then 
become clear that edge corrections to the periodic orbit sum cannot be neglected in this case. Indeed these are 



the only corrections for the particular model that we study numerically in Section IV. Though the absence of edge 



corrections severely affects the spectral determinant, which has a natural cut-off in time, they cancel out in the 



Berry- Tabor formula as discussed in Section III 



II. THE SEMICLASSICAL TRACE 

Relation (|l|) for the eigenvalues becomes exact in the case of the harmonic oscillator actions (||). Likewise, the 
propagator (0), extended to all n, would give the full evolution operator for any given function of the projection 
operators However, the initial semiclassical approximation is to consider the action-angle variables {1,0) as 

power 



appropriate conjugate variables for quantization. From the action representation of the operator [/^, the I 
of Un in (^, the matrix elements in the angle representation turn out to be 



th 



1 

2^ 



N+ 

E 



exp 



9') -I 



(13) 



where we used {e\n) = e^'("+^) V\/2^- 

The Poisson transformation is now applied, changing to the continuous action variable /, so that E is interpolated 
by H{I) as in (0). This results in 



2Trn ^-^ 



d/ exp < — 

n 



I{0 - 0') + 2TTm - ItH{I) 



(14) 



which is exactly equivalent to dl3). Tracing over the angle variables gives 



Tr(C/; 



N) 



E (-1) 



dl exp 



2TrmI - ItH{I) 



(15) 



The exponential in the integrand of (15^ oscillates rapidly in the semiclassical limit, and the largest contributions to 

1/2 — 

the trace, of order h ' , come from the regions near the stationary phase points at 



It 



dHjl) 
dl 



27rm = 



(16) 



which is the same condition as the one given by (||) for the periodic orbits of the classical stroboscopic map. The 
stationary phase evaluation of the trace is therefore tied to the periodic orbits with actions Imj. 

The periodic orbits can now be ordered by increasing period, /. Equation ( |l^ ) also sets the range of possible 
repetitions (M;_, M;+) of the period / orbit. The frequencies, ui = dH{I)/dI in the limited range of action (I_,X+) 
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are also bounded. Thus, the number of periodic orbits in this range increases linearly with /. The full stationary 
phase evaluation of ([l5|) becomes 



Tr([/^)= (-1) 



1/2 



27r 



Mt 



exp^- 



27TmL 



m,l 



(17) 



Not only are the stationary points of the integration specified by the action variables of the periodic orbits, but the 
phase of each contribution is given by the full action in units of h, 



S„ 



It 



dt [pq - H] 



(18) 



except for the geometric "Maslov term", 7r/4. 

Typically, the contributions to the trace from the endpoints J± to the integrals in (|l5|) are only of order h. However, 
note that they cannot be separated from the stationary phase term corresponding to a periodic orbit which lies very 
close to the boundary. In fact, by varying the parameters of the system and /, essentially all the stationary phase 
points pass close to the boundaries. The edge corrections can be obtained as a uniform approximation but this 
erases the simplicity of the periodic orbit expression for the trace. It is shown in the following section that the 
limitation of the allowed phase space to the ring (X_,I+) in no way hinders the Berry- Tabor equivalence. 



III. THE BERRY- TABOR EQUIVALENCE 



A version of the Gutzwiller trace formula that is appropriate to an integrable map is simply obtained by Fourier 
analyzing the traces of in the discrete time 



y e''^Tr((7' ) = Tr = V - 

;=o n 



1 



exp[i(0 — ItE,t) 



(19) 



Here, the eigenvalues of the map are considered as poles of the resolvent (1 — zU^)^^ instead of zeros of det(l — zUn)- 
Though the infinite series of traces in ( p^ ) appears to be a cumbersome alternative to evaluating the finite deter- 
minant, Poisson transformation to the periodic orbit expression can be performed for the semiclassical resolvent. 



Tr- 



1 



1/2 



1 - e^^C/AT 



X exp 



277 



Mt 



hl9 + 2limlm,l — lTH{Im,l) — Ti- 



(20) 



Interpolating the discrete time I again by the continuous time s gives 

/ \ 1/2 



Tr- 



1 



1 - e^^C/iv 



27r 



V 



ds exp 



d^g(J) 



d72 



hs6 + 2ll{mlm,s + flks) — STH{Im,s) — fi- 



The stationary phase points, Sm.kif^), of the integrals in Eq. (21) are given by 

dS„ 



ds 



+ h{9 + 2-Kk) = 



(21) 



(22) 



where Sm,s is the full action of the orbit singled out by condition (|T^). It is important to note that, while in (p^), 
the stationary phase condition defined a periodic orbit of integer period / for the r-stroboscopic map, now I has been 
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replaced by the continuous variable s. This condition can be reinterpreted in two ways: (i) orbits are picked that are 
no longer closed, or (ii) the selected orbit is still periodic, but for a different T{d). In either case, 

^=t^^tE^.s, (23) 
ds dt 

where the energy of the selected orbit appears in the last equality. Thus, the resulting phase of each integral is given 
by the reduced action for the orbits of energy h{6 + 27rfc)/r. The stationary evaluation of the resolvent is 



k.m 



(24) 



where I{E) is the local inverse of H{I). Thus, if h/r is chosen such that the variation of TH{I)/ti does not exceed 2tt 
for / in (X_,X+), there are at most a few orbits contributing to the semiclassical resolvent for each value of 9. Since 
the repetitions must still be summed, the interpretation (ii) for ( p2[ ) as a periodic orbit of another stroboscopic map 
is perhaps the more appealing. 

The condition for the resolvent to be singular is that all the m repetitions are exactly in phase, 

l(ho + 27Tk)] = (j-l)7T , (25) 



which is just the initial Bohr-Sommerfeld quantization. Thus, the periodic orbit evaluation of the map resolvent 
retrieves the full Berry- Tabor equivalence without any need to consider boundary corrections to the trace. The 
reason why this works is that the singularities of the resolvent only manifest themselves by summing the periodic 
contributions for very long times, or multiple repetitions. For these contributions there is an effective increase of the 
large parameter of the stationary phase evaluation of Tr(C/jy) given by (p^. The integrals are then dominated by a 
very narrow region near each periodic orbit, so that the boundary can be ignored for the high repetition of an orbit 
even if it lies almost at the edge. 

The evaluation of the spectral determinant may be regarded as a resummation of the semiclassical sum for the 
poles of the resolvent, such that the (discrete) time of the contributing periodic orbits is cut-off by the dimension of 
the finite Hilbert space. We show in the next section that the error in the traces of the propagator due to the edge 
corrections can be magnified by the spectral determinant. 



IV. MODEL HAMILTONIAN 



To illustrate the calculation of the spectral determinant, we choose the simplest Hamiltonian that is a nonlinear 
function of the action, namely 

H{I) = ^{I-I)' (/>0). (26) 



The periodic orbit condition (m then reduces to 



and the stationary phase evaluation of the trace assumes the explicit form 



27r ^ 



MU'n)^(^) /^5](-l)"cxp 



I / 7rm\ 
-2.™(l+-) 



(28) 



For the case where the parameter T = ti{n+ 1/2) (i.e. one of the quantized actions), Tr([/^) reduces to a "curlicue" . 
These recursively spiraling patterns in the complex plane were shown by Berry and Goldberg to have quite diverse 
characteristics depending on t/H, in the limit N oo. This may also be the case for non-quantized X. 

Two essentially different regions of the eigenspectrum can be studied. In case (i), X > and X± ~ 1 ± ^J2E < 21. 
The energy levels there originate on two different branches of H{I), so there may be near degeneracies. This is a 
similar situation to that expected for an integrable system with more degrees of freedom (see e.g. |l^). Otherwise, in 
case (ii), both X+ and X_ > 21, so that a single branch of the H{I) curve is sampled. The energy levels are then quite 
regularly spaced with approximate separation of Planck's constant times the average dH/dl for this energy interval. 
By choosing E± to satisfy condition (|^) , we also guarantee the regularity of the spectrum of Un itself. 
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FIG. 1: Comparison between the exact (dotted lines) and the semiclassical values (solid lines) of Re[Tr(C/jv)] as a function of 
tjtYi for (a) unweighted and (b) cosine weighted trace. The behaviour is similar for Im[Tr((7jv)] and therefore not shown here. 



Let us first examine the semiclassical approximations to Tr([/^r) for a single period. Notice that the stationary 
phase approximation to the integrals in ( p^ ) would be exact for the quadratic model if the limits could be extended 
to ±oo. In other words, the edge corrections are the only deviation of the periodic orbit approximation with respect 
to the exact trace. These corrections are certainly not negligible for stroboscopic parameters such that a periodic 
orbit is close to one of the edges. The simplest asymptotic edge corrections, as deduced in are presented in 
The computed error due to the edges is displayed in Figure 0. Here we let the time t vary and study Tr(?7Ar) as a 
function ofi/^H, where tn = 2tt/S.N/E is the Heisenberg time. As the time is increased the integrand at the l.h.s. of 
Equation (|l^) has more stationary phase points. The sudden jump in the semiclassical trace in Figure |^ is due to the 
following. For small values of t/tn there is just a single stationary phase point with m — Q. Going beyond t/^H ~ 0.25 
all I TO I < 1 become stationary phase points. (It is near such transition points that the semiclassical approximation 
is worse.) This situation repeats itself regularly with increasing i/^n adding more stationary points to the sum ([2^). 
As seen in Figure |l]a, the semiclassical trace as given by ( p8| ) is not an accurate approximation to the exact Tr(C/7v)- 
The suppression of the edge correction by means of a cosine weighted trace, namely 



A. Case (i) 




(29) 
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FIG. 2: Comparison between the exact (white histograms] 
and (b) Re[Tr(C/j^)] as a function of / for = 48. 



and the semiclassical (grey histograms) values of (a) Im[Tr(C/^)] 



improves the agreement enormously, as shown in Figure 0b. Our calculations have been performed for case h = 
0.01, E = 30.4571694, and I ^ 20.3489573 giving TV = 1 521. 

We now return to the original programme of obtaining the eigenvalues of Un, defined for a fixed t, by resumming 
periodic orbits. The first step is to calculate the traces of C/^. The coefficient of the characteristic polynomial are then 
easily computed by (|ll|). The roots of the resulting polynomial are obtained numerically. The precision in finding the 
roots, limits this method to polynomials with TV < 50. We chose the model Hamiltonian parameters accordingly. 

To illustrate case (i), we take the parameters E and I to be the same as above and h = 0.35, which gives TV — 48. 
Figure |^ displays the real and imaginary parts of Tr(J7^) for ? < TV using the semiclassical expression ( p8| ) contrasted 
with their exact values. As expected, the semiclassical approximation works reasonably well for the lowest values of I 
and gives very poor results for the largest ones. It is worth noticing that similar calculations with the addition of an 
imaginary part to the time in the trace ( p7| ) show a dramatic improvement in the agreement between the semiclassical 
and the exact traces. This is further evidence that the differences are due to edge corrections and not to numerical 
inaccuracies. 

Next, we calculate the coefficients Ck of the characteristic polynomial (^ using the recursion relation (|lTl). The 
symmetry relation ( p^ ) allows us to use only traces corresponding to times shorter than Tii/2. In Figure ^ we compare 
the exact coefficients Ck with those resulting from the semiclassical traces, with and without symmetrization. To 
better illustrate the whole range of Ck we show a situation where TV = 28, corresponding to h = 1.0, E = 100.1234, 
and T — 41.2345. As before, the semiclassical approximation works well for the low k coefficients {k < 5) and. 
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FIG. 3: Comparison between the exact (gray), the semiclassical (white) and the symmetrized semiclassical (black) values of 
(a) Im(cfc) and (b) Re(cfc) as a function of k for A'^ = 28. 



obviously, also for the corresponding higher fc's if we impose the symmetry (g^. For the remaining coefficients, the 
agreement with the exact Cfc is poor. 

Finally we arrive at the eigenvalues of C/^. Figure || shows the results for % = 0.35, = 30.4571694, and 2 — 
20.3489573, which is a rather typical situation. The exact roots of the characteristic polynomial lie on the unit circle, 
as they should. The standard semiclassical approximation destroys unitarity and the roots of the corresponding 
characteristic polynomial are no longer restricted to the unit circle. The enforcement of the symmetry ( p^ makes 
I det(?7)| = 1. As a consequence the roots either lie exactly on the unit circle, or appear in pairs, one inside and one 
outside the circle. More precisely, the self-inversive symmetry of the characteristic polynomial renders a symmetry 
in its zeros: if Zk is a root then l/z^ is also a root. Indeed, it has been shown by Bogomolny et al.||l^ that on 
average only about 57% of the roots of a random symmetrical polynomial lie on the unit circle. Unfortunately, upon 
symmetrization individual roots do not necessarily come closer to the exact ones, as compared with the standard 
procedure. Figure 4 shows that they can even be pushed out, unless the standard semiclassical roots are already close 
to the exact ones, in which case symmetrization tends to improve the accuracy. Unfortunately this behaviour is not 
very systematic. 
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FIG. 4: Comparison between the exact (asterisks), the semiclassical (open circles) and the symmetrized semiclassical (open 
squares) roots Zi of the characteristic polynomial for A'^ = 48. 



B. Case (ii) 




We switch now to case (ii) and consider the situation where the energy E > or equivalently, both X+ and 

X_ > 2T. As mentioned before, in this situation the levels are almost equally spaced with separation of TidH/dl. 
From the semiclassical point of view, the most important distinction to case (i) is the absence of the m = orbit. 
Actually, for our simple model there will be no tori with orbits of period / < N/2 in case (ii). Indeed, the fraction, /, 
of zero traces in the semiclassical approximation is easily seen to be: 



(30) 



The optimum of / w N/2 is approached for E- w 0, which implies X < 0, so we choose fi = 1.0,1 = —1.03489573 
with E_ = 0.265678 and E+ = 1250.0, which corresponds to TV = 48. Keeping r as defined by (0) gives Tr(C/^) = 
for I < 24. The results are summarized by Figure ||. 

As in case (i) the simple semiclassical approximation shows roots outside the unit circle. The standard semiclassical 
improvement by symmetrization ("bottom up symmetrization" ) erases all system information, since there are only 
nonzero traces for long times. Indeed, the spectral determinant ( [l0| ) reduces to 

Pn{z) = 1 + i-lf dctiUN)zN. (31) 

Instead, one can maximize the semiclassical information by symmetrizing the lower coefficients in Pn{z) from the 
higher c^'s , containing the traces with long orbits ("top down symmetrization"). However, the final result is no 
better than the one without symmetrizing at all. 

By shifting the parameters in case (ii) so that H{I) becomes almost linear, within the interval (X_,X+), the 
agreement between the different approaches becomes much more reasonable than in Figure ^ However, if one recalls 
that the mean level spacing is known and the levels are almost equally spaced, some caution is then required. The 
exact traces lead to a characteristic polynomial with self-inversive symmetry. For small values of /, where there is 
no stationary phase and the semiclassical traces are zero, the (modulus of) exact traces are small compared to unit. 
On the other hand, we observe that, for large / values, the semiclassical approximation fails to reproduce the exact 
traces with the same precision. Thus, the best semiclassical agreement corresponds surprisingly to the case where all 
traces are taken as zero (bottom up). As shown in Fig ||, the top down symmetrization is superior to the standard 
procedure only where the eigenvalues are sufhciently accurate. It then brings most of the roots back to the unit circle. 
However, since it employs some inaccurate large I traces, (top down) symmetrization also produces pairs of roots lying 
at opposite sides of the unit circle, like in case (i). 
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FIG. 5: Comparison between the exact (crosses), the semiclassical (open circles), the "bottom" symmetrized semiclassical (open 
squares) and the "top" symmetrized semiclassical (open triangles) roots Zi of the characteristic polynomial for N = 48. Panel 
(a) displays all roots, while (b) is a scale blow up of the same data showing that the nice agreement is only apparent. 



V. CONCLUSIONS 



The main assets of using spectral determinants to resum the periodic orbits in the Gutzwiller trace formula for 
chaotic quantum maps are: (a) the requirement of unitarity is partially incorporated, (b) the necessity of handling 
very long (and exponentially many) periodic orbits with periods larger than T///2 is eliminated and (c) the method 
produces isolated levels rather than a smoothed density of states. Nevertheless, the method remains problematic in 
general, because of the necessity of accounting for a formidable number of periodic orbits in the semiclassical limit, 
which is a daunting task from the classical point of view. 

Even though semiclassical spectra of integrable systems are obtained most efficiently by generalizations of the 
Bohr-Sommerfeld rules, one could expect that spectral determinants would still be more successful tools for integrable 
spectra than the trace formula, since they take some account of unitarity and here the classical periodic orbit structure 
is much simpler than for chaotic systems. Our study shows that this is not the case. 

We find large errors in the semiclassical traces, which could, in principle, be completely fixed by accounting for 
edge corrections. Indeed, this has been verified in the case of the first order corrections in the Appendix. However, 
to follow such a path would be at odds with the spirit of the present work. Rather than solving a trivial model, our 
purpose was to assess the efficacy of the spectral determinant, symmetrized or not, as a tool of reducing the effect of 
inaccuracies of the semiclassical traces in the calculation of eigenvalues. The spectrum of our model is exactly known 
and it allows complete control over each stage of the approximation. Thus, we have worked solely within a framework 
that is based on the usual amplitudes of periodic orbits in the trace formula. Even though our map has been obtained 
by looking stroboscopically at a simple system with continuous time and its boundaries are entirely arbitrary, it is 
arguable that our results may resemble those for a quantum map that results by taking a "Bogomolny section" of a 
Hamiltonian system j^]. The corresponding classical map must also have a boundary and it will also be integrable, 
if this is a property of the original Hamiltonian system. 

The errors in the semiclassical traces, discussed in Section IV, yield inaccurate coefficients for the characteristic 
polynomial, rendering poor approximations for the spectra of integrable systems. The symmetry property of the 
spectral determinant does not necessarily lead to better approximations. Actually in most of the cases studied there 
is only improvement where the unsymmetrized eigenvalues were already reasonably accurate. We can now understand 
this to account for the encouraging results obtained by Creagh |12| for the perturbed cat map. Even though the exact 
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eigenvalues are not determined semiclassically as in our case, his system lies very close to the linear map, where the 
periodic orbit traces are exact, so that the symmetrization always fixes the eigenvalues on the unit circle. In contrast, 
if one perturbs the traces in the present integrable system continuously from their exact values to their semiclassical 
approximation, the eigenvalues may collide on the unit circle and be knocked off as a pair, of which we see many 
examples in our model. We have also tried fitting the spectrum by treating det Un as a free parameter, without any 
essential qualitative improvement of the results. 

It is important to mention that the edge correction may be neglected if the eigenvalues are obtained by Fourier 
transformation over long times, so that the Berry- Tabor equivalence does not need to be "dressed" because of the 
boundary. Of course, this is not much help if the system is not integrable. Indeed, an important point concerning 
the Bogomolny approach is that no assumption is made about characteristics of the dynamics. Recently, this method 
has also been extended successfully to describe the eigenstates themselves of a chaotic system . The present study 
of the integrable limit, though not obtained by a surface of section, suggests that caution may be required in any 
attempt to extrapolate the general section-method to nearly integrable, or mixed systems. 
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APPENDIX A: EDGE CORRECTIONS 



The asymptotic form for the edge corrections follows by basic application of the expansion 

0(|z|-3) . 



,2 e i 



Ciz)+iS{z)^ / dte'^' - ^sign(z) 

Jo V2 TTZ 



(Al) 



where the first term gives the stationary phase contribution; actually, it is the complex conjugate with z real which 
turns out to be needed. The second term is responsible for the edge corrections. Expanding the argument of the 
exponential in (O) to quadratic order in /' = / — ; for the Hamiltonian in (|2q), gives 



Arg — — 2t:vilp , 

Ti \ It 



27r^m^ It _^i2 

^ y 



(A2) 



Rescaling the action variable by y'uh/lT matches the argument of (Al) with that of ( A2). After a little algebra, the 
edge corrections e take the form 



2'n:'^rn? 



1 ( .iT^ 
— exp —t — 1, 

1+ ^\ 2h + 



1 f It ^2 



It ^-^ 1 ?i \ It 

For a discussion of edge corrections for the general case where the phase in the integral (hSl) is not quadratic, see 



(A3) 
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